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Proton Stability in Six Dimensions 
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We show that Lorentz and gauge invariance explain the long proton lifetime within the standard 
model in six dimensions. The baryon-number violating operators have mass dimension 15 or higher. 
Upon TeV-scale compactification of the two universal extra dimensions on a square T 2 /Z2 orbifold, 
a discrete subgroup of the 6-dimensional Lorentz group continues to forbid dangerous operators. 



The stability of the proton is one of the most intrigu- 
ing problems in physics. The current limit on the proton 
lifetime is 1.6 x 10 33 years if the dominant decay mode 
is p — > e + 7r°, and 1.6 x 10 25 years independent of the 
decay mode |j]] . The conserved quantity associated with 
proton stability, the baryon number B, is the charge cor- 
responding to a global U(1)b symmetry, or one of its 
discrete subgroups. It is important to explain why B 
is conserved to such a high degree of accuracy. In this 
letter we demonstrate that the spatial symmetry of two 
universal extra dimensions, accessible to all the standard 
model fields, naturally constrains the proton lifetime to 
be greater than the experimental limits. 

The minimal standard model, including operators of 
dimension four and lower, automatically conserves B. 
This theory, however, is valid only up to some scale 
Asm where new physics becomes manifest. The effect 
of physics at scales above Asm is described most gen- 
erally by higher dimension operators, suppressed by the 
appropriate power of Asm- Among these are B violating 
operators of dimension six: (q c L qi / )(q c L lr,), (d R uii)(q c L li,), 
(Q C LQL){u c R e R ), (u c R d R )(e R u R ), where only one Lorentz 
invariant form is displayed for each four-fermion term; 
the weak-eigenstate quarks, qL, Ur, d R , and leptons, It,, 
e R , may belong to any generation. The experimental 
limits on the proton lifetime impose a bound on the coef- 
ficients of these operators: CsM/Ag M < 10~ 24 TcV~ 2 . If 
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SM 



> 



10 16 GeV, then the coefficients Csm can safely be 
as large as order one at that scale. However, the natural- 
ness of the Higgs sector in the standard model restricts 
Asm to lie in the TeV range. The Csm's must therefore 
be extremely small — below 10 -23 or so. 

Various explanations for the smallness of B violation 
have been proposed. In the minimal supersymmetric 
standard model, a discrete symmetry is needed to banish 
i3-violating terms of dimension three and four from the 
superpotential. The discrete symmetry may be a rem- 
nant of a 17(1) gauge group 0], or else one could hope 
that it is preserved by quantum gravity. In the case of 
large extra dimensions, the smallness of Csm has been 
suggested to arise from the localization of the quarks and 
leptons at different points inside a thick brane |l . 

Here we concentrate on the chiral 6-dimensional stan- 
dard model, in which all standard model fields propa- 
gate in two extra dimensions. The current bound on the 
compactification scale of two universal extra dimensions 



is 1/R ^, 500 GeV 0], suggesting a rich phenomenology 
fej. The standard model in two universal extra dimen- 
sions is especially appealing because of the properties 
of the Lorentz group in six dimensions, as well as the 
constraints imposed by anomaly cancellation. The 6- 
dimensional standard model is an effective theory valid 
up to a scale M s « 5/R, above which the 6-dimensional 
SU(3)c x SU(2)w x U(1)y gauge interactions become 
non-perturbative. 

In six dimensions, the 50(1, 5) Lorentz symmetry has 
two irreducible spin-1/2 representations generated by 
S Q ' 3 /2, where S Q/3 = i[T a ,T' 3 ]/2 in terms of the an- 
ticommuting 8x8 matrices T a , a = 0,1,. ..,5. The 
product r°r 1 ...r 5 has eigenvalues ±1 defining the two 6- 
dimcnsional chiralities. The irreducible gauge [SU(3)c x 
SU(2)w x J7(l)y] and gravitational anomalies cancel for 
only two chirality assignments, up to an overall sign 



U- , V. 




(1) 



where Q+ , U- , X>_ and £± , E^ , A/ip are the 6-dimensional 
quarks and leptons, respectively, and a generational in- 
dex is implicit. We assume that only one of these two 
assignments applies to all fermion generations, since then 
the SU(2)w global anomaly demands 3 mod 3 genera- 
tions [pi. Each fermion generation includes a gauge sin- 
glet field, Afzf, such that the gravitational anomaly can- 
cels. The 6-dimensional chiral quark and lepton fields 
have four components, and if two dimensions are ap- 
propriately compactified on an orbifold, then their zero- 
modes may be identified with the left- and right-handed 
standard model fermions QL,UR,dR,lr,,eR, as well as 
three right-handed neutrinos vr. We first describe the 
constraints on B violation from 6-dimcnsional Lorentz in- 
variance, and then show that the constraints are equally 
tight even with the lesser symmetry remaining after orb- 
ifold compactification. 

We begin by observing that SO(l, 5) has an 50(1, 3) x 
£7(1)45 subgroup: the combination of the 4-dimensional 
Lorentz symmetry associated with the x^ coordinates, 
jU = 0, 1, 2, 3, and the symmetry under rotations of the 
x 4 , a; 5 coordinates. Any chiral 6-dimensional fermion, 
S&±, may be decomposed under this subgroup as 
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TABLE I. Transformation properties of the quark and lep- 
ton fields under SO(l, 3) x [7(1)45, and the corresponding 
zero-modes under orbifold compactification. 



where L and R are £0(1,3) chiralities, defined by the 
projection operators 



P±l = Ptr = \(1T^) 



(3) 



Since S™' 3 /2 are the generators of the spin-1/2 repre- 
sentations of SO(l, 5), the 50(1, 3)-chiral fermions have 
11(1)45 charges given by the eigenvalues of E 45 /2: +1/2 
for *± L , and ±1/2 for V± R . In Table I we list the 
charges of the quarks and leptons. The covariant deriva- 
tive, D a , transforms as a vector under 50(1,5), so that 
D4 + iD 5 has 1/(1)45 charge — 1. 

To construct the _B-violating operators in six dimen- 
sions, we first note that since all quark fields carry 
B = +1/3, SU(3)c gauge invariance demands that for 
any such operator the number of quark fields (minus 
the number of charge-conjugated fields) is a multiple of 
three. Only operators constructed of fields that have zero 
modes can induce at tree-level the B-violating processes 
searched for in experiments so far. From Table I we see 
that all quark fields containing zero modes have £7(1)45 
charge —1/2, so that an operator with \AB\ > 1 has 
a 11(1)45 charge given by (— 3/2)AB plus the sum of 
the charges of all lepton fields included in that operator. 
Hence, all B -violating nucleon- decay operators allowed by 
Lorentz invariance in six dimensions must involve at least 
three quarks and three leptons. 

Using this key fact, it is straightforward to find the 
lowest dimension operators, invariant under 50(1, 5) and 
standard-model gauge transformations, that (after com- 
pactification) can induce B- violating nucleon decays. We 
consider first the £+ chirality assignment. The oper- 
ators then first appear at mass-dimension 16 in the 6- 
dimensional theory, each involving six fermions and one 
covariant derivative p = DpT°: 

(C+V-) 2 (Jf-pV-) , (S-0V-) (Jf-T a V-) 2 , 
(£ + V_)(7J-Q+)(ir-pV_) , 
(AT_ J pW_)(A7_r Q P_) 2 , (JJ_Q + f (JJ_pV_) , (4) 



where we have exhibited only the Lorentz covariant bi- 
linears with the smallest number of Y matrices; permu- 
tations of p and T" are also allowed. The 50(1,3) chi- 
ralities of the zero modes imply that only the terms in 
(0) which do not contain T 4,5 can induce nucleon decay. 
Each of the above dimension-16 operators enters with a 
coefficient proportional to M~ 10 leading to a strong (and 
adequate) suppression of B violation after compactifica- 
tion to four dimensions. We estimate their effects after 
first considering dimension-17 operators, whose contribu- 
tions to proton decay turn out to dominate. 

There is only one proton-decay operator of dimension 
17 (modulo Fierz transformations, and the insertion of 
two or three S"' 3 matrices), invariant under 50(1,5) 
and standard model gauge transformations (in particular 
U(1)y), that does not involve the gauge singlets J\f~: 



017 = 



° 17 'Z+V-fH, 



M. 



(5) 



where H is the charge-conjugated Higgs doublet in six 
dimensions. Dimension-17 operators involving 7V_ have 
a similar form, and will be discussed below. All other 
-B-violating nucleon-decay operators require more deriva- 
tives, 7i fields, or fermion bilinears. 

Upon integrating over the compact dimensions, x A , a; 5 , 
restricting the fields to their zero modes, and replacing 
the zero-mode Higgs doublet by its VEV (vh ~ 174 GeV), 
On gives rise in the 4-dimensional theory to the operator 
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where A45 is the area spanned by x and x 5 . For the 
T 2 /Z2 square orbifold of radius R constructed in [H, 
A45 = 2ir 2 R 2 . This operator is non- vanishing only if the 
generational indices of the three £ + 2?_ bilinears in On 
are not all identical. Hence, On induces proton decays 
into e~ ir + tt + vv or \i~ ir + ir + vv with the two neutrinos be- 
longing to different generations. Final states with more 
pions or a K + are kinematically suppressed. The proton 
width is given by: 
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(^) n ^F(nn), (7) 



where m p is the proton mass. If the process occurs pre- 
dominantly via the decay of the d valence quark [p|, 
whose constituent mass m p /3 is used to set the scale in 
Eq. (R), then the form factor F(tttt), which describes the 
probability for the two pions to be formed, is expected 
to be of order unity. The kinematical phase-space factor, 
$5, defined here as the dimensionless quantity that re- 
lates the squared amplitude to the decay width, is tiny for 
a five-body decay. We estimate $5 % (47r)~ 7 x O(10~ 4 ), 
where the factors of 47r account for angular integration 
and the additional suppression is due to integration over 



the magnitudes of final state momenta. The process 
p — > /!~7r + 7r + z^ has an even smaller $5. Thus, Oyj 
yields a finite proton lifetime: 
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12 
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(8) 



where the quantities in square brackets are of order one 
or larger. There are no published experimental limits for 
five-body proton decays. For comparison, the searches 
for p — > e~7r + 7r+ || set a limit of t p > 3 x 10 31 yr. 
Thus, for C17 < O(l), the proton lifetime is orders of 
magnitude longer than the current experimental limits. 

The On operator also induces i3-violating neutron de- 
cays, including n — > e~ r K + vv^ [i^ix^vv. The four-body 
decays proceed via the fusion of two d valence quarks, 
which is suppressed by their wave function overlap in- 
side the neutron ||. This is compensated though by the 
larger constituent mass, ~ 2m p /3, of the dd pair. Com- 
pared with T(p — > e~Ti + ir + vv) given in Eq. ([7|), the neu- 
tron partial width in these modes is enhanced by $4/ $5, 
where $4 < (47r) -5 x O(10~ 3 ) is the four-body phase- 
space factor. Thus, the inverse width of the 7?-violating 
neutron decays induced by On is smaller than t p by three 
orders of magnitude. The experimental limit on the par- 
tial mean life of n — > e~7r + is 6.5 x 10 31 yr M]. If a 
dedicated search for n — > e~Ti + vv were to yield a com- 
parable limit, then only a weak bound of Gyj %, 0(1) 
would be imposed, with uncertainties especially due to 
the sensitive dependence on M s and R. 

These are striking results. B- violating nucleon decays 
are adequately suppressed even with the scale of B viola- 
tion in the TeV range, providing only that the coefficient 
C17 is not larger than order unity. In fact, Cn is likely to 
be significantly smaller than unity because the operator 
017 is composed of chirality flipping bilinears and there- 
fore can be expected to arise with strength proportional 
to small Yukawa couplings. 

We now return to the other operators. For the 
dimension-16 operators of Eq. (|J), the presence of the 
covariant derivative leads to an additional factor of m p /3 

— 1/2 

in the 4-dimensional amplitude in place of a A 45 Vh/M s 
factor. Thus, there is a supression factor of 10 3 in the 
widths relative to those arising from dimension- 17 opera- 
tors. As we mentioned above, all dimension-17 operators 
other than 0i7 involve at least one A/1 field. The tiny 
neutrino masses suggest that the couplings of its zero 
mode, v a, are small. Thus, the coefficients of these op- 
erators may naturally be substantially smaller than Cyj. 
Nonetheless, some could be relevant because they allow 
new decay modes: p — > ir+vvv, n — > ir°vi>i>, and other 
final states with mesons and three neutrinos, [n — » vvv 
requires a three-quark fusion S whose small probabil- 
ity is not compensated by the gain in phase space.] The 
experimental limits on such modes are of order 10 32 yr 



A/1 field, (£ + D_) 2 (A7_Q + )H, and a few with more A/1 
fields. The bounds on the coefficients of these operators 
are as loose as those on C17. 

Consider next the £_ chirality assignment. It leads to 
a different set of proton-decay operators, invariant under 
S0(l,5) and gauge transformations: the dominant ones 
are six-fermion operators of dimension 15, 

1 r -1 /_„ \ 2 



M? 



C 15 (£ + U.) + C' 15 (C_Q, 



C_T a U- 



(9) 



Other six-fermion operators involve A/+ fields, but again, 
their coefficients may be small. The operators (pi) lead 
to proton decay into three anti-leptons (one or two are 
positively charged) and a number of mesons. For the 
dominant decay mode, p — ► e + e + 7r~j?, we estimate 
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where F(yt) is a form factor of order unity. Despite the 
reduced phase space, the process p — > /i + e + 7r - £/ is more 
constraining due to better data on decays into muons 
§: T p /Br{p -> ii+X) > 10 31 yr, where Br(p -► /j+X) 
is the branching fraction for inclusive decays. We de- 
rive a constraint, C\§ % 10~ 2 (i? x 0.5 TeV) -5 , which is 
rather loose given that C15 may naturally be as small as 
some Yukawa couplings. Nonetheless, the prospects for 
observing B violation look better compared with the C+ 
chirality assignment. Recall that no known theoretical 
argument determines which of the two chirality assign- 
ments of Eq. (TO) is preferred. 

We have so far considered proton decay arising from 
operators that respect the SO(l, 5) symmetry. But since 
the compactification of two dimensions breaks 5*0(1, 5), 
including its 1/(1)45 subgroup, we must next study 
whether the reduced symmetry allows operators that 
could induce proton decay at an unacceptable level. A 
simple and symmetric choice for the compactification of 
and a; 5 dimensions is a T 2 /Z% orbifold of equal 



the x 

radii || . The 
27r_R in the x 4 
tions, ifiix 4 , x 5 



'torus" T is taken to be a square of size 
x 5 plane, with periodic boundary condi- 
= ip(x 4 + 2irR, x 5 ) = <p(x 4 , x 5 + 2irR) for 
any field ip. The square has a Z A symmetry, being invari- 
ant under ir/2 rotations in the x 4 ,x 5 plane. Now, since 
the generator of the £7(1)45 rotations acting on fermions 
is E45/2 [the spin-1/2 representation is "double-valued" 
under 27r rotations], upon compactification on a square, 
the [7(1)45 symmetry is broken down to a Zg group whose 
elements are given by U n with U = cxp[i(7r/2)S]45/2] and 
n = 0,l,...,7. 

Compactification on T 2 would allow only vectorlikc 
zero-mode fermions in the 4-dimensional theory, so it is 
necessary to introduce the observed SO(l, 3) chirality of 
the standard model fermions by orbifolding the square. 
Each field is taken to be either even or odd under the Z2 
orbifold transformation, defined by a n rotation in the 
There exists an operator of this type with a single x 4 ,x 5 plane. The assigned Z 2 parity of the quarks and 




FIG. I. The square T 2 /Z 2 orbifold of radius R: the funda- 
mental region is shaded, and the four fixed points under 
Z2 are marked by •. 

leptons can be read from Table I: a field has a zero mode 
if and only if it is even. The fundamental region in the 
x 4 , x 5 plane of the square orbifold is shown in Fig. I. 

The important point is that the square T 2 /Z% orb- 
ifold does not break further the Z§ symmetry of the 4- 
dimensional effective theory. This is because the Z% orb- 
ifold projection commutes with the Z% transformation. 
A fermion ij) of charge z under Z% transforms as follows: 



Uip(x 4 ,x 5 



k iz-K jl 



a e 



iP(\x 5 



(11) 



where k = 1 (fc = 2) if ij) is odd (even) under Z^. The 
transformation of ip(x , x 5 ) depends on the sign of x 5 , 
through a = +1 (a = -1) for a; 5 < (x 5 > 0). Note 
also that the presence of orbifold fixed points (see Fig. 1) 
does not break the Zg symmetry. 

Having established that the residual symmetry of the 
effective theory is SO(l, 3) x Zg, where Z$ is the sub- 
group of {7(1)45 defined above, it follows from SU(3)c 
gauge invariance and from the charges given in Table I 
that all operators constructed only of zero-mode fields 
satisfy the selection rule §Af? ± 5 AX = mod 4. In 
particular, no proton-decay operators with less than 6 
fcrmions are allowed and our previous operator analy- 
sis remains valid (with the dominant B- violating effects 
given by dimcnsion-17 (-15) operators in the £+ (■£-) 
chirality assignment). Furthermore, it is also clear that 
in this model there are no n—n oscillations (which require 
a AB = 2 operator) and no Majorana masses (AX = 2). 

The remarkable suppression of B violation established 
here is due to the Zg symmetry. But if, for example, 
x A and x 5 were compactified on a rectangular orbifold 
(T 2 jZi with different radii), then the Z% would be re- 
duced to a Z±. It would follow that dimension- 10 op- 
erators such as (Q + T a Q + )(Q + T a C + ), which do not in- 
clude terms with only zero-modes, would induce proton 
decay via loops. [For the £_ chirality assignment, the 
leading X?-violating operators have dimension 12, e.g. 
(Q C + T a Q + ){C c _T a V^)H.} It is therefore important to 
explore whether the Zg symmetry is indeed left intact by 
the compactification. Here we offer only a few comments 



regarding the choice of vacuum. It has been noted that 
the square configuration is an extremum of the Casimir 
energy of bulk fields, including gravity, while the rectan- 
gular one is not jllj . Although the lowest order compu- 
tation indicates that the square configuration is a saddle 
point of the effective potential (the minimum corresponds 
to a "rhombus" that has a Z§ symmetry, still forbidding 
X?-violating operators with less than three quarks and 
three leptons), it could be that a complete computation 
will reveal this configuration to be a minimum. Alterna- 
tively, it is possible to freeze the shape of the torus to be 
square by a Z± orbifold identification, in which case the 
Z^ symmetry remains unbroken. Either of these possi- 
bilities would leave intact our conclusions based on the 
T2/Z2 square orbifold regarding B violation. 

We finally mention that some of the fields describing 
gravitational fluctuations are also charged under the Zg 
(the components of the metric along x 4 and a; 5 ). How- 
ever, we need not worry about the effects of these fields 
for nucleon decay if the radius stabilization mechanism 
makes them heavier than the proton or if they are pro- 
jected out by a Z± orbifold identification. 

In conclusion, we have shown that the combination of 
standard-model gauge invariance, 6-dimensional Lorentz 
invariance, and compactification of the two extra dimen- 
sions on a T 2 /Z2 orbifold of equal radii naturally sup- 
presses proton decay to acceptable levels, even with the 
scale of baryon number violation in the TeV range. 
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